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The stationary gravitational field of two identical counter-moving beams of pure radiation is 
found in full generality. The solution depends on an arbitrary function and a parameter which sets 
the scale of the energy density. Some of its properties are studied. Previous particular solutions are 
derived as subcases. 
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CN ■ I. INTRODUCTION 

O ; 

' The gravitational field of incoherent light is induced by the energy-momentum tensor T^^ = /ii k^ki^ of pure radiation 
(null dust) with energy density fii and travelling along the null vector k^. The study of steady beams and pulses of 
^ pointed radiation has a long history. In the linear approximation to general relativity this was done about 70 years ago 
^\ [0,11 ■ Later, Bonnor gave exact solutions for steady and time-dependent beams |g[ which belong to the well-known 
class of pp- waves (algebraically special solutions of Petrov type N and no expansion) [Q. He showed that beams 
00 ■ shining in the same direction do not interact. He also gave solutions with charged null dust |^ and spinning null dust 
1^ . Counter- moving light beams interact non- linearly, some explanation being proposed in Ref . . Particular exact 
solutions were found only recently Their Killing vectors necessarily have twist to prevent the focusing of rays. 

One of these solutions has been also obtained as a limiting case of the double-dust solution Later, the general 
double-dust solution was described flli . 



> 



X 



[ In this paper we find the stationary gravitational field of two identical colliding beams of light in full generality. 
■ Some of its properties are studied and the previous results are rederived as special cases. 
0^ ' In section 2 the properties of the energy-momentum tensor are studied, the metric and the system of Einstein 
. equations are written down when the energy densities of the beams coincide. Some simplifications are done next and 
' the main equation is obtained. In section 3 a new radial variable is introduced and two methods for the derivation 
— ..^ , of the general solution are given. The explicit formulas depend on an arbitrary function and a free parameter. The 
^ consequences of elementary flatness are studied. In section 4 an algorithm for simple algebraic solutions is presented 
y [ and the Kramer solutions are derived as an example. The difficulties in the second method, which contains the 
5h 1 important case of constant energy density, are clarified. In section 5 some of the properties of the general and the 
particular solutions are studied. The expansion, the acceleration and the vorticity of k'^ are given in the most general 
J> metric. A criterion for the appearance of closed timelike curves (CTC) is given. The static solution of Kramer does 
not have CTC but its stationary counterpart does have them. Section 6 contains some discussion. A comparison is 
made between the general double-dust solution and the general solution for colliding light. 



II. ENERGY-MOMENTUM TENSOR AND FIELD EQUATIONS 

We shall work in the stationary formulation of the problem when the cylindrically symmetric metric is given by 

SI 

ds^ = -e^" {dt + Adifif + [e^'^ {dx^ + dz^) + W^dip^] (1) 

and there is only radial cc-dependence. The static metric may be obtained by complex substitution; t — > iz,z — s- 
if, A —t I A. The energy-momentum tensor of the two beams reads 

Tf_,i, = iJ.ik^_,ky + H2ltj.lu, (2) 

fc^fc^ = 0, = 0, (3) 
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k'^h. = -1. (4) 

Here are the energy density profiles of the beams and fc'', are their 4- velocities, which are null and normalized 
as in Eq. (4). As will be shown later, there is only one solution with the minimal set of two non-zero components. 
On the other hand, naturally fc^ = P = 0. Therefore, let us introduce 

fc* = = V, k'f = 1'^ = $, k' ^-V- ^ Z (5) 

and T = V -\- A'^. The beams move in the positive or negative z-direction. For stationary metrics they should have 
the same twist otherwise the two conditions in Eq. (3) lead to contradiction. In the static picture the twists are 
opposite. The obvious relation = and Eq. (4) give 

g2«j,2 _ g-2«yf^2^2 ^ g2(fc-«)^2^ (6) 
^2urp2 _ g-2up^2^2 ^ g2(fe-«) ^2 ^ ^ ^j-^ 

Their combinations express Z and $ in terms of the other quantities 

Z^^e^-\ (8) 



^2^2 ^ g2u |^g2uy2 _ _J _ ^g^ 

The covariant conservation of the energy-momentum tensor (which follows from the field equations) reads 

where g is the metric's determinant. In our case ~ Q and the first term vanishes. The second yields 

(e2«) V2 ^ 2 (^e2"y _ (g-2u^2 _ g2u^2y $2 ^ 3 (fc - u)' e^C^-")^^, (11) 

where ' means x-derivative. The restriction $ = immediately gives the particular solution k = 2m, as announced 
above. 

To write the Einstein equations we use the combinations of Ricci tensor components utilized in Rcfs. |11 l^ , [l3| , [l^ . 



We shall work, however, with curved mixed components R'^. The only non-trivial ones are the diagonal components 
and i?^, Ri . The Einstein equations for the corresponding Tj^ components give after some rearrangement 

^ = 87rpe2(^-"), (12) 
2^ - + ^ = 0, (13) 

ii'W p4:U A/2 

4tt A '2 

k" + u'^ + -tS^ = 16^pe2/c-4«p^2^2 (^5) 
AW ^ ' 



W 



iiirpe^^W'^T. (16) 



Here 2p = + ii2- Units are used with G = c — \. The equations for the rest components are satisfied either 
identically or when (/ii — (12) Z = 0. Therefore, we accept in the following that pi — ^2 = P- There are 7 equations 
(8,9,12-16) for 8 unknowns u, k, W, A, p, V, $, Z. Eq. (13) gives an expression for A' 
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A' ^±2e-^^wVB, B = u'^-^, (17) 

which allows to eliminate A' from Eqs. (14-16). The difference between Eqs. (14,15) leads to a relation between k — u 
and W 



1 - W 



e 



W ' W 



(18) 



An unimportant constant has been set to 1 in the first of these equations by rescaling the x-coordinate. Now Eqs. 
(12,18) show that p can be expressed through W 

Snp^WWe"^! (19) 

or W can be found as a solution of a rather complicated equation if p was given 

{WW")' - 2W" - (In p)' WW" = 0. (20) 

Eq. (8) defines Z in terms of W, Eqs. (14,9) express T and <& in terms of u and W. Thus all characteristics of the 
problem have expressions as functions of or m and W. The arbitrary function, however, can be only one and a 
relation between them follows from the crucial Eq. (16). While Eqs. (14,15) are not sensitive to the sign of A' in 
Eq.(17), Eq.(16) should be solved for both signs of A' , leading to different solutions. They are encompassed by taking 
the squares of both sides of Eq. (16) and inserting T^, <i>^ and A' from Eqs. (14,15,17) to obtain 

This equation involves only u and W due to Eqs. (17,18). From the definition of T it follows that if A, <i>, T is a 
solution of Eq. (16), so is —A, — T. In other words, the equations are invariant to a change in the direction of 
dragging and the sign of the twist of the pure radiation's 4- velocity.. This is different from the vacuum case ||l^ and 
the case of perfect fluid [|l3|, where the r.h.s. of Eq. (16) vanishes and the direction of rotation is not significant. In 
our case the two signs of A (or <i>) lead to a double valuedness of the energy density p, as will be seen from Eqs. (29, 
37) below. 

III. THE GENERAL SOLUTION 

Let us introduce the new radial variable p = W' and denote the p-derivatives with a subscript. Eq. (12) shows that 
this is possible when p ^ 0- Let us introduce also the important function X 

X = WW" = (22) 
Another fundamental quantity is /, introduced by 



Wp 



Then B becomes 



S--^, b = f-pf+p'-p, (24) 

where 6 > in order to have a real A' . This constrains to some extent the arbitrary function /. After the use of the 
last three equations the main Eq. (21) turns into a quadratic equation for X 

CiX^ + C2X + 0-3 = 0, (25) 
Ci=bl- 4bfp ifp - 1) , (26) 
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C2=4b[{2f-p)bp-2b{2fp-l)], 



(27) 



{2f-pf-4b 



Every function / corresponds to two solutions 

-C2±VD 



2Ci 



£» = C| - 4C1C3. 



The coefficients Cj are simplified when the definition of b is used 

Ci=p(4-3p)/2 + 2(3p/-2/-p)/p + (2p-l-/f 



= 46 [p (4 - 3p) /p + 3p/ - 2/ - p] , 



C3 = 462p (4 - 3p) 



(28) 

(29) 

(30) 
(31) 

(32) 



It is interesting that only 4 combinations of p and / enter these expressions. Another surprising fact is that D factors 
out into terms that do not contain fp 



D = 166^ 



(3p/ -2f-pf-p{A- 3p) (2p - 1 - ff = 6463 {p - 1) (3p - 1) . 



(33) 



X must be real, hence D > 0. This restricts the p coordinate outside the interval (1/3, 1). The position of the axis 
is at p = 1, because W {x) ^ x when a; ^ 0. Thus the realistic interval is given by p > 1. It is also obvious that Eq. 
(25) is quadratic with respect to fp 

p{A- 3p) X^f^ + [2 {3pf -2f-p)X^+ 4p (4 - 3p) bX] fp + {2p-l- ff X^ 
+4 (3p/ ~2f~p)bX + Ap (4 - 3p) 6^ 0. (34) 

After these remarks we are ready to present the general solution of the problem of two colliding light beams. There 
are two methods to do this. The easier one is to choose an arbitrary function / {p) and find the two X from the root 
formula (29). Then Eqs. (23,18,22,17,19,8,9 and 14) yield respectively u,k,W,A,p,Z,^ and V 



u = j -^dp, k = j ^ "^^ ^ dp + fco. 



W = Woe 



P Vb 



We-'^'^dp, 



v2 



2 _ 2b + fpX . 

2X ' 

$2 ^ 2b+ {fp - l)X ^2u-2j y = T - A$. 

2Wq X 

The radial metric component Qxx becomes gpp 



9pp — e 



2{k-u) 



(35) 
(36) 
(37) 
(38) 
(39) 

(40) 



One can further pass from p to the physical length L, so that Qll = ^ 
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L = J Va^dp ^ Woe! . (41) 

The constant Wq is positive by definition. The hniits of the integral for A and the integration constant fco are 
determined by the elementary flatness condition 

lim ^ (e~2"M/2 _ ^2u^2^ 1/2 ^ (42) 

W and A should vanish at the axis. One can see from Eq. (17) that A = o (W) when x ^ even if there is a 
pole in B. Therefore, it can be neglected in Eq. (42) like in many other problems pl] , p^ , p^ . The function / (and 
b) also should vanish at the axis to keep u' regular. Using the I'Hopital rule and passing to p, Eq. (42) becomes 
(2p - / - 1) e~'' 1 when p ^ 1. This gives A: (1) = which fixes fco. Inserting W and p from Eqs. (36,37) into Eq. 
(12) supplies the relation 



SirpoW^ ^ 1. (43) 

It indicates that po is positive. Eq. (43) has been used in Eqs. (38,39). Eqs. (37,43) show that X > 0, which 
places another restriction on /. Furthermore, the r.h.s. of Eqs. (38,39) should be positive, leading to the inequality 
26+ (/p — 1) X > Q. Finally, we have kept A positive in Eq. (36), allowing for a change in the sign of $ instead. It is 
determined by the sign of the l.h.s. of Eq. (16), equivalent to the sign of (4/ — 2p) b + bpX. 

The second, more natural but more difficult method, is to provide the sum of the density profiles 2p. Then from 
Eq. (20) or (37) it follows that X satisfies a linear first-order differential equation 

Xp-(lnp)pX-2 = 0. (44) 

Its boundary condition states that X (1) — 0, due to W (1) = 0. Its solution is 

X = 2pnX (45) 
Ji P 

The next step is to determine / from Eq. (34). This is, in fact, a first-order, highly non-linear equation, whose 
solutions seem to exist only in numerical form. Its boundary condition is / (1) = 0. After X and / are found, the 
rest of the task is fulfilled by Eqs. (35-40) where the two methods coincide. 

The general solution thus depends on one arbitrary function / (satisfying several positivity requirements) or a 
positive p and the free parameter po- The situation is similar to the general double-dust solution |ll| . The constant 
Pq sets the scale of the energy density. The form of the solution is determined by Eq. (25) or its equivalent Eq. (34). 
They do not contain or any other parameters. An explicit solution needs the explicit calculation of only three 

integrals, /f, /if^ and j ij^. 

IV. PARTICULAR SOLUTIONS 

Let us search for some simple concrete solutions. Following the first method, let us choose / in such a way that D 
becomes a square leading to an algebraic X. According to Eq. (33) we must have 

b={p-l)Cip-l)b\, (46) 

where 61 is a polynomial. Eq. (24) shows that the simplest possibility is constant bi and linear /. There are four 
such cases 

h = l-p. k = 2u, (47) 

/2 = ^, k = -u, (48) 

/3 = ^, k = 3u-2liiW, (49) 
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/4 = 2p-l, k = \nW. (50) 

Then Eq. (29) yields 8 solutions for X 

Y I ^^ Y P(p-l)(3p-l)(4-3p) 

^ia = (p-l)(3p-l), Xu = 3p2 _ 4p + 2 ' ^ ^ 

_ (p-l)(3p--l)(4-3p) _ p(p-l)(3p-l) 

g^^^ , ^2.- 2^^^ , (52) 

-'^Sa = —X2b, Xsb = —X2a, (53) 

Xia — ~Xih, X4h — —Xia- (54) 

The first four coincide with solutions I-IV from Ref. |^]. As was mentioned already, realistic solutions have positive 
X, at least in some region around the axis. Setting p = + 1 we see that Xia and X2b are positive everywhere, while 
Xib and X2a are positive for 1 < p < 4/3. These solutions are suitable for interior solutions. The other four solutions 
are physically unrealistic. X41, and X^a are negative everywhere, while X^a and are negative on the axis and in 
its vicinity. This fact was noticed in Ref. ||^. 

We have shown that when ^ = 0, k — 2u. The opposite is not true since Xn, leads to non-vanishing Thus the 
solution with /i and Xi is the simplest possible solution. It was thoroughly investigated It is worth to point 

out how it arises from the general formalism. When k = 2u, u also picks up an integration constant fco/2, while the 
one in p becomes —k^. Eqs. (35-40) give 

u^Y-^'\ ^"^ = ¥-^1^ e^'^^2"\ (55) 
= W^R^Y-'^/^, A^2^/^WoR^, (56) 
p = 2^'^poY^, $ = 0, (57) 

Z ^T^ -^Y^/^, (58) 
v2 



where 



y = i(3p-l) = l + ^i?2. (59) 



2 ' " ' 2 

These results coincide with the results of Refs. [^,^ when the identifications 

22/3 [2 

87rpo = — A^, ^=Y3^^ (60) 

are made. The free parameter A reflects the free scale of the energy density. 

Eq. (46) sets a combinatorial problem from polynomial algebra, which probably has additional solutions with 
more complex /. Furthermore, when / and X are given analytically, the integrals in the metric may be calculated 
numerically and a solution, partially explicit, partially numerical, obtained. 

The second method for generating solutions also has an analytical start by choosing p. Then X is found compara- 
tively easy from Eq. (45). Here are four examples 

pi — const, Xi = 2 (p — 1) , (61) 

P2=p-", X2 = ^(p-p-"), (62) 
a + 1 

P3=c-p, = 2(c-p)ln^, (63) 

c — p 

P4 = e-^ X4 = 2 (1 - e^-P) . (64) 

The constants fulfil c > l,a > 0. The first case is the simplest from a physical point of view - collision of beams 
of constant energy density up to some po, where a matching to the vacuum stationary solution should be done. 
Unfortunately, even for such simple X Eq. (34) is so complicated that we have been unable to solve it even numerically. 
This hopefully can be done on a more powerful computer platform. 



6 



V. SOME PROPERTIES OF THE GENERAL SOLUTION 



We have calculated with the help of GRTensor the characteristics of the vector field fc'^ in the general stationary 
cylindrically symmetric metric (1). The shear's expression is too lengthy to be displayed. The expansion O, the 
acceleration and the vorticity vector read 

e ^W~'^ [{W<i>)' ~2W^{u' - k')] , (65) 
a' ^<f(v' + 2u'V+ ^ y (66) 



(67) 



M/2 

+ {u - k') {Z^ - $2) + (68) 
a" = $ [2 {k' -u')Z + Z'\ , (69) 

2u 

= — [2 (2m' ~ k') AVZ + ZAV' + ZA'V - Z'AV] , (70) 

lu'P = -^[2{2u' -k')VZ + V'Z -VZ'], (71) 
^ 0, (72) 

When (f> = these formulas coincide with those from Ref. Both and are proportional to <f> and vanish in 
this case except for a^, which becomes 

It also vanishes for the general double-dust solution and for the only case fc = 2u of colliding light with $ = 0. The 
motion of the source is therefore geodesic. When $ 7^ Eq. (8) shows that 

= (fc' - u') (75) 

which vanishes only in the vacuum case k ^ u. The expansion does not vanish automatically when the general solution 
is inserted in Eq. (65). Vorticity does not depend on $ at all and Eqs. (70-73) coincide with the first relations in 
Eqs. (28-30) from Ref. 

It is well known that CTC exist when g^^ < 0. We have 

g^^ = (e-4"iy2 _ ^2) (7g) 
and the general solution possesses CTC when 

-2mti/ , 2«T 



e-2«VF< / -^e-'^'Wdp. (77) 
'1 X 



For the concrete solution (55-60) this yields 
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which is satisfied when r > \/3/A, so that there are CTC. In the static case Y = 1 — \^r^ and inequahty (78) becomes 

r2r-2/3 + IaV^^ < 0. (79) 

The l.h.s. is always positive and consequently there are no CTC. This is also true for the exterior solution - the Levi 
Civita metric. In the stationary case this role is played by the Lewis solution ||T^. It possesses three distinct classes 
and only the Weyl class is locally, but not globally static. It seems possible to induce the Lewis class of the exterior 
solution when po is big enough. This class contains CTC. This discussion of exterior solutions is valid also for the 
general solution. However, the vacuum solution is not a subcase of the general solution. It has p = 0, ly = a;, so that 
p is not a variable but a constant. 



VI. DISCUSSION 



We have obtained in this paper the general cylindrically symmetric and stationary (no time dependence) metric 
when T^y is given by Eqs. (2-4) and — The components of the null vectors and the energy density are found 
simultaneously with the metric. Although the solution looks quite complicated, minimal assumptions have been made 
along its derivation. Diagonal metric does not allow colliding solutions, hence, the reason for the appearance of A. It 
prevents the focusing of rays and the appearance of singularity [|||^ . In the case of dust two non-zero components of 
are enough for a general solution with arbitrary energy density . In the case of light they give a solution with 
specific density profile and a third component $ is necessary. In the dust case the field equations decouple and can be 
solved one by one, giving, however, an involved expression for the energy density in terms of the arbitrary function. 
In the light case there is an additional quadratic non-linearity, because the square of Eq. (16) must be taken in order 
to encompass all solutions. It reflects further in the quadratic algebraic Eq. (25) for X or the quadratic differential 
Eq. (34) for fp. One can choose arbitrary positive p and flnd X easily. The expression for / and for those parts of the 
metric depending on it, however, can be only numerical. The other option is to choose /, but then there are several 
inequalities to be satisfied for a realistic solution. 

An important property is the traditional appearance of CTC in rotating cylindrically symmetric solutions, which 
forms one of the differences between the stationary and the general static solution. This takes place both in the 
interior and the exterior. 

In conclusion, this paper can be considered as a further step in the study of the collision of the gravitational fields of 
beams of coherent or incoherent light. The award motivating the pursuit of this issue is the existence of singularity-free 
colliding solutions ||^. 
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